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The cholesterol-dependent cytolysins (CDCs) attack cells by punching large holes in their 
membranes. Lectinolysin from Streptococcus mitis is unique among CDCs due to the pres- 
ence of an N-terminal lectin domain that enhances the pore-forming activity of the toxin. 
We recently determined the crystal structures of the lectin domain in complex with various 
glycans. These structures revealed the molecular basis for the Lewis antigen specificity of 
the toxin. Based on this information we have used in silico molecular modeling to design 
a mutant toxin, which we predicted would increase its specificity for Lewis y, an antigen 
found on the surface of cancer cells. Surprisingly, we found by surface plasmon resonance 
binding experiments that the resultant mutant lectin domain exhibited higher specificity for 
Lewis b antigens instead. We then undertook comparative crystallographic and molecular 
dynamics simulation studies of the wild-type and mutant lectin domains to understand the 
molecular basis for the disparity between the theoretical and experimental results. The 
crystallographic results revealed that the net number of interactions between Lewis y and 
wild-type versus mutant was unchanged whereas there was a loss of a hydrogen bond 
between mutant and Lewis b compared to wild-type. In contrast, the molecular dynamics 
studies revealed that the Lewis b antigen spent more time in the binding pocket of the 
mutant compared to wild-type and the reverse was true for Lewis y. The results of these 
simulation studies are consistent with the conclusions drawn from the surface plasmon 
resonance studies. This work is part of a program to engineer lectinolysin so that it will 
target and kill specific cells in human diseases. 

Keywords: cholesterol-dependent cytolysins, Lewis antigens, molecular dynamics simulations, protein engineer- 
ing, surface plasmon resonance. X-ray crystallography 



INTRODUCTION 

Lectinolysin (LLY) is a pore-forming toxin derived from some 
strains of Streptococcus mitis and S. pseudopneumoniae and is 
a member of the cholesterol-dependent cytolysin (CDC) family 
(Farrand etal, 2008). These pore-forming bacterial toxins are 
secreted as soluble monomers that assemble on the surface of 
cholesterol- rich cell membranes. The resultant pores are com- 
posed in excess of 30 monomers, and are greater than 150 A in 
diameter (Hotze and Tweten, 2012). Perforation of the membrane 
by these pores results in cell lysis. The sequence similarity is high 
(40-80%) between CDC family members. 

Streptococcus mitis can cause diseases such as infective endo- 
carditis and septicemia (Hall and Baddour, 2002; Huang etal, 
2002; Gowda etal, 2003; Kennedy etal, 2004). Serum iso- 
lates from Kawasaki disease patients were found to contain S. 
mitis human platelet aggregation factor (PAF or Sm-hPAF), so 
named because it was found to affect light scattering properties 
of human platelets, thought to be as a consequence of platelet 
aggregation (Ohkuni etal., 1997). Farrand etal. (2008) found 



that the Sm-hPAF gene sequence encodes a predicted CDC struc- 
ture (Rossjohn etal., 1997; Polekhina etal., 2005), most closely 
related to intermedilysin (Farrand etal., 2008). Having estab- 
lished that PAF was in fact a member of the CDC family, it was 
renamed lectinolysin, because of its ability to bind carbohydrate 
(see below). 

While LLY shares a number of characteristics typical of CDCs, it 
has a unique 162 amino acid N-terminal domain (LLY^^^; Farrand 
etal., 2008). Sequence comparisons show that the LLY^^*^ domain 
shares significant identity with fucose-binding proteins. Glycan 
array experiments revealed that LLY was highly selective for the 
difucosylated glycans Lewis y (Le^) and Lewis b (Le^) antigens 
(Farrand etal., 2008). Small-angle X-ray scattering (SAXS) mod- 
eling predicts that LLY^^^ is joined to CDC domain 1 so that the 
Le antigen-binding site is exposed on the outer surface of the pore 
(Feil etal, 2012). 

Lewis antigens are blood group determinants with very rigid 
structures (Yuriev et al, 2005). Le^ is a type 1 antigen and adsorbed 
onto the surface of blood cells. Le^ is a type 2 antigen which, in 
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healthy individuals, is expressed at low levels in tissues includ- 
ing epithelial cells and serum (Yuriev et al, 2005; Westwood et al, 
2009). However, expression is elevated on the surface of a wide 
range of epithelial tumor cells including colon, lung, and ovary 
(Sakamoto etal, 1986; Miyake etal, 1992; Yin etal, 1996). Thus, 
Le^ is a distinctive tumor marker and therefore a promising target 
for directing drugs to cancer cells. 

Various approaches have been taken to produce Le^ -specific 
cancer therapeutics, including the development of Le^ -specific 
humanized monoclonal antibodies (Kitamura etal., 1994; Clarke 
etal, 2000a,b; Scott etal, 2000), and anti Le^' T cells (Westwood 
etal., 2005). A chimeric approach has also been taken whereby a 
Le^ -specific monoclonal antibody was conjugated to the CDC, lis- 
teriolysin O, to form an immunotoxin (Bergelt et al., 2009). LLY is a 
naturally occurring, potential oncotoxin, with in-built tumor tar- 
geting specificity for the Le^ antigen. However, cross reactivity with 
Le^ antigen decreases its usefulness as a specific, cancer-targeting 
molecule. 

Previously we have determined the crystal structure of LLY^^*^, 
and its complexes with fucose and Le^ and Le^ antigens (Feil et al, 
2012). The overall fold of the LLY^^^ domain exhibited similarities 
with fucolectin domains of the CBM family 98 glycoside hydrolase 
of Streptococcus pneumonia and Anguilla anguilla agglutinin. The 
binding of Le^ and Le^ to LLY^^^ was very similar with the excep- 
tion of an additional hydrogen bond between Le^ and amino acid 
residue Tyr62 (Feil etal., 2012). 

Using these structures, we have explored the possibility of 
mutating LLY^^^ to develop a Le^ -specific binding domain. In the 
work presented here a Y62H mutation was designed and con- 
structed with the aim of removing the hydrogen bond between 
the AT-acetylglucosamine (GlcATAc) of Le^ and LLY^^^ and instead, 
gain a hydrogen bond with Le^. However, binding studies 
described here show that LLY^^*^Y62H has a greater affinity for 
Le^ over Le^, whereas LLY^^^wt has a preference for Le^ over 
Le^. In order to explore the molecular basis for this unexpected 
finding we determined the crystal structure of LLY^^*^Y62H in 
complex with Le^ and Le^. We then compared these structures 
with those of the wild-type lectin domain complexed with the 
same Lewis antigens. The crystallographic studies were comple- 
mented by molecular dynamics simulations in order to explore 
the importance of protein motion and interaction with the ligands 
over time. 

MATERIALS AND METHODS 
INSILICO DESIGN OF MUTANTS 

The crystal structures of the lectin domain of LLY in complex 
with Le^ antigen (PDB code: 3LEK) and Le^ antigen (PDB code: 
3LEG) were visually inspected using Pymol^. In silico mutations 
were created and assessed also using the Pymol program. Figures 
were created in Pymol and using Chemaxon software^. 

SURFACE PLASMON RESONANCE BIOSENSOR BINDING ANALYSIS 

All surface plasmon resonance (SPR) experiments were performed 
at 25° C using a Biacore T200 instrument (GE Healthcare). Direct 
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assays were performed with Le^ and Le'^ antigens (Sigma) injected 
over mutant and wild-type LLY^^^ domain proteins immobilized 
on a CM5 chip (GE Healthcare). Immobilizations were performed 
in 1 X HBS-P running buffer [10 mM HEPES, 150 mM NaCl, 
0.005% (w/v) Tween 20]. Binding experiments were performed in 
1 X HBS-P + running buffer [10 mM HEPES, 150 mM NaCl, 
0.05% (w/v) Tween 20] containing 1 mM MgCl2, 1 mM CaCl2, 
and 0.2 mg/ml bovine serum albumin. 

LLY^^^wt and LLY1^^Y62H proteins were immobilized in two 
separate channels on a CM5 chip using a standard amine cou- 
pling protocol. Briefly, the chip surface was activated with 
a single 5 min injection of freshly prepared 1:1 50 mM N- 
hydroxysuccinimide:200 mM 3-(iV,iV-dimethylamino)propyl-N- 
ethylcarbodiimide. Protein coupling was achieved by three 5 min 
injections of LLY^^Vt or LLY^^^Y62H solution (50 |xg/ml in 10 mM 
sodium acetate, pH 5.0). To deactivate residual reactive sites, 
lectin domain coupling was followed by a 5 min injection of 1 
M ethanolamine (pH 8.5). Approximately 3000 response units 
(RU; 1 RU = 1 pg of protein/mm^) of LLY^^Vt (channel 2) and 
4000 RU of LLY1^^Y62H (channel 3) were coupled. Channel 1 was 
activated and blocked, as above, for use as a reference surface. 

Le^ and Le^ were injected, in duplicate, in a twofold dilution 
series from 1 mM to 7.5 |xM over the immobilized wild-type and 
mutant lectin domains. Thirty-second injections of Le antigen 
at 40 |xl/min were followed by a 60-s dissociation period. Base- 
line returned to 0 almost immediately after the association phase 
ended, so regeneration of the surface was not required. 

Binding data were processed and analyzed using Scrubber 2 
software^ (version 2.0c). Rapid association and dissociation rates 
made data fitting to a kinetic model and subsequent calculation 
of kinetic rate constants and impractical. Consequently, 
equilibrium dissociation constants (Kd) were derived by fitting 
binding responses at equilibrium to a 1:1 steady- state affinity 
model available within Scrubber. 

MUTAGENESIS, EXPRESSION, AND PURIFICATION OF THE LLY'^'' 
DOMAINS 

A mutant of LLY^^S LLY1^^Q190C, with an N-terminal 6 x His 
tag and a TEV protease cleavage site for His-tag removal, was 
used in this study. The LLY^^*^Q190C coding region corresponds to 
LLY residues 38-190 (GenBank accession number AB051299.1). 
The C-terminal residue was mutated to cysteine (Q190C) for 
protein labeling studies. The point mutant displays wild-type 
activity and is located far from the Lewis antigen-binding site 
(Farrand et al, 2008). Henceforth, LLY^^^Q190C will be described 
as LLY^^^^wt in this work. Expression was carried out in E. coli 
BL21, and the protein purified using nickel resin and size exclu- 
sion chromatography as previously described (Feil etal., 2012). 
LLY^^Vt was mutated with QuikChange™ Site-Directed Muta- 
genesis Kit (Stratagene) to create the Lewis antigen-binding site 
mutant LLY1^^Y62H. LLY^^^wt and LLY1^^Y62H were purified to 
>95% purity, as determined by SDS gel electrophoresis. Pro- 
teins were concentrated to 10 mg/ml in 10 mM Tris-HCl pH 
7.2, 10 mM NaCl, and 5 mM dithiothreitol, and were stored 
at -80°C. 
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CRYSTALLIZATION OF THE LLY'^'^YGZH MUTANT 

All crystallizations were performed using the hanging drop vapor 
diffusion method at 21°C. 2 |xl of protein was mixed with equal 
volume of precipitant and hung over 0.5 ml of well solution. 
LLY^^^Y62H crystallization conditions were determined by fine 
screening around the conditions that proved successful for the 
wild-type protein (Feil etal., 2012): 2 M MgS04 and 100 mM 
Tris-HCl buffer pH 8.2-9.0. Optimization included the addition 
of the Hampton additive screen (Hampton Research, CA, USA) 
10% (v/v) in the hanging drops. The optimal crystallization condi- 
tions were 2.4 M MgS04 and 100 mM KCl with 100 mM Tris-HCl 
pH 8.75-9.0. Le^ and Le^ antigens were soaked into LLY^^^Y62H 
crystals by adding 1 |xl 1 M MgS04, 100 mM Tris-HCl pH 9.0, 
100 mM KCl, 5 mM Le^' or Le^ to the crystal drop for 1 h at 21°C. 
Soaks were carried out immediately prior to cryoprotection and 
flash freezing in liquid nitrogen. Crystals were cryoprotected for 
X-ray data collection by adding glycerol in increments of 5% (v/v) 
to a final concentration of 20% (v/v) to the crystal drops. 

crystallographic studies of the lly'^'^ygzh mutant 

Diffraction data were collected at the MX2 beamline at the Aus- 
tralian Synchrotron in Clayton, Victoria. The data collection 
was controlled using Blue-Ice software (McPhillips etal., 2002). 
The diffraction data were processed using the HKL2000 suite 
(Otwinowski and Minor, 1997). Model building was performed 
with COOT (Emsley and Cowtan, 2004) using the published 
wild-type crystal structure of LLY^^*^ (PDB code: ILEO) as a 
starting model. Data were refined with REFMAC 5 (Murshu- 
dov etal., 1997) from the CCP4 program suite (CCP4, 1994). 
Restrained positional and isotropic temperature factor refinement 
was employed until convergence and the refinement was mon- 
itored using the i^free residual. The Lewis antigens were fully 
occupied in the binding sites of both structures and their tem- 
perature factors were very similar to those of the surrounding 
side-chains. 

MOLECULAR MODELING OF THE LLY'®'' INTERACTION 
WITH THE Le'^ AND Le^ ANTIGENS 

To investigate the mechanism of association and disassociation 
between LLY^^^ and Lewis antigens four molecular models were 
used: two being the published crystal structures of LLY^^^^wt bound 
to either Le^ or Le^ antigen (Feil etal., 2012) and the other two 
being the crystal structures of LLY^^^Y62H, complexed to the same 
ligands, that are described here. The molecular dynamics program 
NAMD (Phillips etal., 2005) was used to model ligand disasso- 
ciation of each complex with a 30-ns simulation from the initial 
starting bound conformation, repeated 100 times for a cumulative 
simulation of 3 |xs for each system. Models were initially solvated 
with TIP3 water under periodic boundary conditions of dimen- 
sions 48 A X 48 A X 64 A. Charges were neutralized with NaCl 
for a total ionic concentration of 150 mM. The bound calcium 
ion seen in the crystal structures was included in the simulations. 
Dynamic molecular modeling was conducted at a theoretical pH 
of 7.4. Each simulation run started with an equilibration phase 
of 0.5 ns, where the protein backbone and Lewis antigen was har- 
monically constrained to their starting positions at 310 K and 
an NPT ensemble (an isothermal and isobar ic ensemble where 



the number of moles (N), pressure (P), and temperature (T) are 
conserved). Subsequent 30-ns production simulations were run 
with no constraints at 310 K under NVT ensembles [a canonical 
ensemble where the number of moles (N), volume (V), and tem- 
perature (T) are conserved]. Trajectory snapshots were captured 
every 100 ps. At the completion of 100 repeat simulations, all pro- 
duction runs per model were consolidated into a single trajectory 
and the frames root-mean-square (rms) deviation were centered, 
based on the ELY alpha carbon protein backbone. The Lewis anti- 
gen was then subject to clustering analysis, whereby the top 20 
ligand conformational clusters with a cut-off of 1.5 A were gen- 
erated and averaged. Ligand clusters with a rms deviation value 
less than 3.5 A from the initial position were considered to be in a 
bound conformation. 

RESULTS 

IN SILICO DESIGN OF MUTANTS 

The difference between Le^ and Le^ antigens is the core disaccha- 
ride linkage (Gaip 1 -4 versus Gaip 1 -3) with an opposing projection 
of the iV- acetyl glucosamine (GlcNAc) iV- acetyl and -CH2OH 
groups of the iV-acetylglucosamine monosaccharide as shown 
in Figure 1. The overall Lewis determinants maintain the same 
conformation in both free and bound states. 

In the crystal structure of LLY^^*^ wt the GlcNAc of Le^ is within 
hydrogen bonding distance (3.3 A) of Tyr62 (Figure 2 A; Feil 
etal., 2012). This interaction is not able to occur with Le^, thus 
suggesting that mutating Tyr62 may alter the affinity for Le^. 
Various mutations of this residue were created in silico. Muta- 
tion to a histidine residue would maintain structural similarity 
of the binding site whilst also creating the potential for a new 
interaction between the histidine side-chain and Le^ as illustrated 
in Figure 2B. Replacement of Tyr62 with His was expected to 
remove the interaction between Le^ and the hydroxyl group of 
Tyr62 (Figure 2A). 

SPR BIOSENSOR BINDING ANALYSIS 

Purified LLY1^^Y62H yield was 20 mg/1 E. coli BL21 culture in 
Luria-Bertani (LB) broth. This yield was equivalent to that found 
for LLY^^^wt. The sensorgrams from SPR experiments between 
the Lewis antigens (Le^ and Le^) and either wt or LLY^^*^Y62H are 
shown in Figure 3. The estimated Kj) values are listed in Table 1. 
The data indicate low affinity binding interactions of the Lewis 
antigens with the lectin domains. The wild-type lectin domain 
(LLY^^^wt) has more than threefold greater affinity for Le^ over 
Le^. The LLY^^^Y62H mutant has a 1.6-fold greater preference for 
Le^ over Le^. The mutation causes an overall 2. 5 -fold decrease 
in the lectin domain's affinity for Le^ and doubles its affinity 
for Le^ over the wild-type domain. Overall, the data show that 
wild-type LLY^^^ preferentially binds Le^ antigen while the mutant 
LLY1^^Y62H exhibits higher affinity for Le^ than for Le^. 

CRYSTAL STRUCTURES OF LLY'^'^YBZH IN COMPLEX 
WITH Lev AND Le'^ ANTIGENS 

LLY1^^Y62H crystallization was optimized in 100 mM Tris pH 
8.75-9.0, 2.4 M MgS04, and 100 mM KCl. Small bipyrami- 
dal crystals of dimensions 0.15 mm x 0.12 mm x 0.12 mm 
appeared in 6-8 weeks. After 8 months some crystals had grown 
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FIGURE 1 I Structural comparison of Le"^ and Le^ antigens. Lewis 
antigens Le*-* (A) and Le^ (B) are structurally very sinnilar. They differ in the 



core disaccharide linkages and the opposing projections of their A/-acetyl 
(blue arrow) and CH2OH groups of the GlcA/Ac nnoiety. 





Tyr62 



His62 



FIGURE 2 I In silico LLY'e<Y62H mutant design. (A) Le^ i 

was designed to retain sinnilar 
lines denote hydrogen bonds. 



_ « . . is within hydrogen bonding distance (3.3 A) of Tyr62 fronn the LLY'^^wt domain. (B) AY62H mutation 

was designed to retain similar structural properties to tyrosine, with the possibility that it may also result in a hydrogen bond to the Le^ antigen. Yellow dashed 



to 0.25 mm x 0.2 mm x 0.2 mm. In contrast, LLY^^*^ wt crys- 
tals appeared after 5 days in 100 mM Tris pH 8.2-8.4, 2.0 M 
MgS04 with the same morphology as LLY^^*^Y62H, but reach 
the same size within 10 days (Feil etal, 2012). Lewis antigens 
were soaked into crystals immediately prior to cryoprotection and 
freezing. 

The crystal structures of the LLY^^^Y62H in complex with 
the Lewis antigens was performed as described in Materials and 
Methods. The crystal structures have been deposited in the PDB 
with codes: 4GWJ and 4GWI for the Lewis b and Lewis y com- 
plexes, respectively. Data statistics for the structures are listed in 
Table 2. 

LLY^^^ adopts an eight- stranded P- sandwich fold, composed 
of a five-stranded anti-parallel P-sheet on one side and a 
three-stranded anti-parallel P-sheet on the other side. Three short 



a-helices separate P-strands 1 and 2 and 3 and 4. There is one cal- 
cium ion in the structure that is in the same position as observed 
in structurally related fucolectin domains. There is also a metal ion 
that is CO valently bound to His80, as well as five water molecules, all 
arranged in octahedral geometry. In our previous work we iden- 
tified the ion as either a Ni^+ or Mg^+ ion (Feil etal, 2012) and 
here we have chosen the latter possibility. The Le^ and Le^ antigen- 
binding site is in a cleft at one end of the molecule. A more detailed 
description of the structure has been reported elsewhere (Feil et al., 
2012). In the new mutant structures some additional N- terminal 
residues are observed (three extra in the Le^ structure and one 
extra in the Le^ structure), compared to the published wild- type 
structures (see Table 2). 

The structures of the LLY^^*^wt-Le^ (Figure 4A) and 
LLY^^^Y62H-Le^ (Figure 4B) superimpose very closely, with a 



Frontiers in Immunology | Molecular Innate Innnnunity 



November 2012 | Volume 3 | Article 330 | 4 



Lawrence etal. 



Lectinolysin Lewis antigen specificity 



B 







LLY'^^vt_Le> 



























LLY'^^Y62H_Lc> 





























20 30 

Time 



20 30 

Time 




E 




10-5 10-^ 10-3 

Concentration (ulVI) 



FIGURE 3 I SPR measurements of immobilized LLV'^Cwt and LLY'ecY62H. 

Le^ antigen and Le*-" antigen were injected at concentrations of 7.6-1000 |xM. 
(A) LLY'ec^^ vs LeV, (B) LLYlecY62H vs LeV, (C) LLY'^Cwt vs Le^, (D) 



F 




10 "5 10-4 iQ-3 

Concentration (ul^) 



l_l_YlecYQ2H vs Le*^. Responses at equilibrium fitted to simple 1:1 binding 
isotherms for Le^ (black lines) and Le*-" (blue lines) interacting with (E) 
LLV'^^wt and (F) LLY'®^62H mutant. Duplicate binding data sets are shown. 



rms deviation of the alpha carbon atoms of 0.1 A. A third metal 
ion binding site is observed in the mutant structure and we have 
tentatively identified it as a Mg^+ ion binding site due to the very 
high concentrations of this ion in the crystallization buffer. The 
Mg^+ ion is bound between the carboxylate of Asp97 and the side- 
chain carbonyl of Gln54 and four water molecules, all arranged 
in octahedral geometry around the metal ion. The water structure 



is the same as in LLY wt-Le^ apart from a water molecule that 
is positioned between the hydroxyl group of Tyr62 and fucose 
1 (Fuc 1; Figure 4A). This water molecule is shifted by 1.4 A 
toward the NE2 of His62 in the mutant structure (Figure 4B), 
compared to wild-type (Figure 4A), to optimize its interaction 
with the His side-chain. LLY^^*^Y62H forms the same 19 van der 
Waals interactions with Le^ compared to the wild-type protein: 
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Table 1 | Estimated equilibrium dissociation constants (Ko) for 
LLYiecY62H and LLY'^'^wt domains interacting with Le^ and Le'' 
antigens. 





LLY'ec^t Ko (|xM) 


LLY"ecY62H (|xM) 


Ley 


78 ±9 


202 ±31 


Leb 


234 ± 12 


127 ± 14 



15 with Fuc 1, 1 with Gal, 2 with GlcNAc, and 1 with Fuc 2. 
The B- factor of Le^ is 53.5 in the mutant (overall B- factor of 
protein is 29 A^), compared to 38.2 A^ in the wild-type protein 
(overall B-factor of protein complex is 21 A^). Overall, there is 
no net change in the number of potential hydrogen binding and 
van der Waals interactions between wild-type and mutant Le^ 
complexes. 

Superposition of the structures of the LLY^^*^wt-Le^ and 
LLYl^^Y62H-Le^ gives a rms deviation of the alpha carbon atoms 
of 0.1 A. In addition to the Ca^^ and Mg^"^ ion binding sites 
seen in the wild- type structure there is the additional Mg^^ ion 
binding site as observed in the LLY^^^Y62H-Le^ crystal structure. 
In the structure of LLY^^^wt-Le^ there is one hydrogen bond from 
the hydroxyl group of Tyr62 to the iV- acetyl moiety of Le^ and 
a second one via a water molecule to the hydroxyl group of Fuc 
1 (Figure 4C). In the LLY^^^Y62H-Le^ structure (Figure 4D) the 
bond from Tyr62 to the N- acetyl is replaced by a bifurcated hydro- 
gen bond from His62 to the hydroxyl group of the GlcATAc ring 
and the hydroxyl of the N- acetyl moiety off the ring. Overall, 
there was a net loss of a potential hydrogen bond as the water- 
mediated hydrogen bond between Tyr62 and Fuc 1 seen in the 
LLY^^Vt-Le^ structure has been lost in LLY^^^H62-Le^. The water 
molecules between the protein and Le^ are conserved in wild- 
type and mutant structures. LLY^^^Y62H forms the same 23 van 
der Waals interactions with Le^ as does the wild-type protein: 
15 with Fuc 1, 1 with Gal, 2 with GlcNAc, and 4 with Fuc 2. 
The average B-factor of Le^ is 43.8 A^ in LLYl^^Y62H-Le^ (over- 
all B-factor of protein complex is 23 A^) which is very similar to 
36.7 A^ observed in the wild-type complex (overall B-factor of 
protein complex is 18 A^). Thus, overall, there is a loss of one 
potential hydrogen bond between Le^ and mutant compared to 
wild-type. 

Superposition of the crystal structures of LLY^^^Y62H-Le^ and 
LLYl^^Y62H-Le>^ results in a rms deviation on the alpha carbons 
atoms of 0.1 A. Fuc 1 and the Gal moieties of the Le^ and Le^ 
antigens superimpose quite well whereas GlcNAc and Fuc 2 are 
shifted. GlciVAc in Le^ moves 0.6 A toward the protein due 
to a water-mediated hydrogen bond between the oxygen of the 
GICNAC-CH2OH group to the backbone oxygen of His62. Like- 
wise, Fuc 2 is shifted by about 0.6 A because the same interaction 
is not present in the Le^ complex. 

MOLECULAR SIMULATIONS OF LEWIS ANTIGEN INTERACTIONS 

A total of 30,000 trajectory frames (100 x 30 ns) were collected 
for each Lewis antigen in complex with either LLY^^*^wt or the 
LLY^^^Y62H. In the wild-type simulations, the Le^ antigen resided 
within the binding pocket more than the Le^ antigen (Table 3). 
However, when Y62 was mutated to a His this trend was reversed. 



Table 2 | Crystallographic data processing and refinement statistics 
for LLY'^'^wt and LLY'®'Y62H crystal structures in complex with Le^ 
and Le^ antigens. 





LLY"ecY62H-LeY 


LLY"ecY62H-Le'' 


Data collection 


Space group 


D/l 0 0 

r42A Z 


DA 0 0 


Unit cell dimensions (A) 


67.1, 67.1, 99.4 


66.9, 66.9, 99.4 


Wavelength (A) 


0.95 


0.95 


Tennperature (K) 


100 


100 


Maxinnunn resolution (A) 


1.6 


1.6 


No. of observations 


426,761 


404,032 


No. of unique reflections 


30,563 


29,510 


Redundancy 


14.0 


13.7 


Data connpleteness (%) 


99.2 (93.7) 


96.7 (79.3) 


l/ai 


19.6 (71) 


16.9 (5.3) 


Emerge (%)^ 


11.5 (45.8) 


10.9(571) 


Refinement 






Non-hydrogen atonns 






Protein 


1124 


1103 


Water 


146 


150 


Mg2 + 


2 


3 


^ 9-1- 

Ca^+ 


1 


1 


Ligands 


46 


46 


Resolution (A) 


1.6 


1.6 


f^work (%)^ 


16.9 


18.0 


F^free (%) 


19.1 


21.4 


Rnns deviations fronn ideal geonnetry 






Bond lengths (A) 


0.028 


0.028 


Bond angles (0) 


1.4 


1.4 


Bonded B's 


3.9 


4.1 


Mean B (A^) 






Main-chain 


22.8 


270 


Side-chain 


28.2 


31.9 


Water 


36.9 


29.4 


Ligand 


53.5 


43.8 


Residues observed 


41 to 184 


43 to 184 


Residues in nnost favored regions 


877 


88.4 


of the Rannachandran plot (%) 






Residues in the disallowed regions 


0 


0 


of the Rannachandran plot (%) 







The values in parentheses are for the highest resolution bin (approximately 0. 1 A 
width). ^ Rmerge = ^hkl^ii^i -<l>l/l<l>l, where l; is the intensity for the ith mea- 
surement of a symmetry related reflection with indices h,k,l. ^ Rwork. = ^ll^obsh 
l^calcll^l^obsl' where Fqi^^ and F^^gi^. are the observed and calculated structure 
factor amplitudes, respectively 



with Le residing in the binding pocket more than Le^. These 
results are reflective of the SPR binding data, in that LLY^^^ wt has a 
higher affinity for Le^ over Le^, whereas LLY^^*^Y62H has a higher 
affinity for Le^ over for Le^. 
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FIGURE 4 I LLY'®<^wt and LLY'®<V62H domains in complex with Lewis 
antigens. Comparison between (A) the LLY'^^wt-Le^ complex (orange 
carbons, green Tyr62) and (B) the LLY'®^Y62H-Ley complex (purple carbons, 
green His62) shows that the water molecule shifts toward the histidine 
residue in LLY'®^62H-Ley.The water-mediated hydrogen bond to the hydroxyl 



group ofTyr62 seen in the LLY'^^wt-Le"^ complex (cyan carbons, green Tyr62) 

(C) was lost in the LLY'^^wt-Le'-' structure (magenta carbons, green Tyr62) 

(D) as this residue was replaced by His62. Hydrogen bonds are shown in gray 
and water molecules as black spheres. The fit of the ligands to the relevant 
'^Fq-Fq, electron density map (dark green hash), contoured at la, are shown. 



Lewis antigen dissociation rates were plotted as a function 
of occupancy (defined as <3.5 A deviation from the original 
binding site) for all 100 simulations (Figure 5). This indicated 
a rapid dissociation rate which supports the SPR data where 
the on and off rates were too fast to be calculated. There is 



Table 3 | Cumulative totals of bound Lewis antigen over the course 
of 100 independent, 30 ns simulations starting from the Lewis 
bound conformation. 



wt LLY'e^62H 

Ley 17,260 (57.5%) 13,835 (46.1 %) 

Le^ 14,824(49.4%) 16,080(53.6%) 

All binding events witiiin 3.5 A rms deviation of tine original antigen bound 
conformation were considered "bound." 



an overall trend for the antigens to occupy the binding site less 
as the simulations progressed. However, it is important to note 
that the curves are not purely exponential and contain distinct 
areas of peaks and troughs. This is indicative of the ligands mov- 
ing rapidly in and out of the binding site of both LLY^^^wt and 
LLY1^^Y62H. 

The simulation data reiterates the SPR finding that mutation 
of LLY^^^Tyr62 to His62 has changed the binding specificity from 
a preference for Le^ by wild-type to a specificity for Le^ by the 
LLY1^^Y62H mutant. 

DISCUSSION 

The high level of expression of Le^ on the surface of epithelial 
tumor cells and low level expression elsewhere makes cell-bound 
Le^ a highly attractive target for anti-tumor agent delivery. In addi- 
tion, the rigid nature and distinctive structure of Le^ lends itself 
to the ready development of highly specific protein recognition 
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FIGURE 5 I Molecular dynamics modeling of the association of Lewis 
antigens with the LLY'®*^ binding sites. Traces are averaged for all 100 
sinnulations. Overall, the Lewis antigens occupied the binding site of either 
the wild-type or LLY'®^62H less as the sinnulation progressed. However 
there are numerous areas of re-association, resulting in a spike in 
occupancy. 



domains, whether in antibodies or other glycan-binding 
proteins. 

Previously, we pubHshed the crystal structures of LLY^^^^wt in 
complex with the Le^ and Le^ antigens (Feil etal., 2012). These 
structures, together with computational modeling, were used 
to design and construct a glycan-binding domain with altered 
specificity for the Lewis antigens. The LLY^^^Y62H mutation was 
suggested due to its structural similarity of His to Tyr and the 
potential to create a new interaction with Le^. 

Surface plasmon resonance binding studies showed that 
LLY^^^^wt has a higher affinity for the Le^ antigen over Le^. This 
relationship was reversed for the designed LLY^^^Y62H mutant, 
i.e., the Le^ antigen had higher affinity than the Le^ antigen. These 
results were unexpected, as the molecular modeling predicted that 
mutation of Tyr62 to His should decrease the number of hydro- 
gen bonds with Le^ and increase the number of potential hydrogen 
bonds with Le^. The expectation was that the mutation would lead 
to an increase in affinity for Le^ and decrease in affinity for Le^, 
but in practice we found the reverse to be true. 

To experimentally determine the interactions formed in the 
binding site as a result of the LLY^^^Y62H mutation, we solved 
the crystal structures of LLY^^*^Y62H in complex with the Lewis 
antigens. In the LLY^^^Y62H-Le^ complex, there was no addi- 
tional bond formation between the Lewis antigen and His62, only 
a shift of a water molecule toward the His, as compared to its 
position relative to Tyr62 in the wild-type (Figures 4A,B). In the 
Le^ complex there was a net loss of a potential hydrogen bond to 



His62 (Figures 4C,D). In summary, the crystal structures did not 
show the bonding as predicted from the in silico mutant model 
(Figure 2) but revealed a net decrease in hydrogen bonding inter- 
actions to Le^ and no net change in bonding interactions with Le^. 
However, the SPR data revealed that the mutant had an increased 
affinity for Le^ and a decreased affinity for Le^. 

Crystal structures are a static average snapshot of what happens 
in solution, so we used molecular dynamics studies to simulate the 
lectin domain-Lewis antigen binding mechanism in solution. In 
these simulations it was observed that the amino acid residues 
around the Lewis antigen binding site are mobile, particularly 
Tyr62. Thus, the inferred hydrogen bonding interactions between 
protein and ligands observed in the crystal structures may not 
necessarily be persistent in solution. This may help explain the 
unexpected consequences of our mutation. 

The molecular dynamics studies provided us with comple- 
mentary results to the SPR data: both approaches showed that 
the association and dissociation rates with Lewis antigens are 
extremely rapid, and that the relative affinities between antigens 
and wild-type or LLY^^*^Y62H domains are generally weak. The 
SPR data showed the strongest interaction being that of LLY^^^^wt 
to Le^ with a Kj) of 78 |xM and the weakest between LLY^^^^wt and 
Le^ with a Kd of 234.4 |xM. Relative affinities of the wild-type or 
LLY^^*^Y62H binding sites for antigens was measured in the sim- 
ulations by calculating how much time the Lewis antigens spent 
in the binding pocket. These dynamics studies agreed with the 
general trend of the biological results, suggesting that there was 
a reversal of Lewis antigen affinity when Tyr62 was mutated to a 
His. Interestingly, and again in concurrence with SPR data, the 
dynamics studies also suggest that the Lewis antigens bind weakly 
and transiently to the LLY^^*^ domain. 

A recent example of protein engineering to alter Lewis antigen 
specificity is that of Noro virus virion protein 1 (NoV VPl), which 
specifically binds Le^ (Kubota et al., 2012). However, NoV VPl has 
a deeper binding site with more extensive contacts to the Lewis 
antigen structure than seen in LLY^^^. A single residue (Gin) in the 
Le^ binding site was mutated (to Asn) with the aim of increasing 
the affinity of the protein for Le^. The mutation increased the 
width of the binding site, allowing bonding between both fucose 
rings and amino acid side-chains. As expected, the mutant protein 
displayed a higher affinity for Le^. 

In contrast, designing a high affinity Le^ -specific LLY^^^ mutant 
is more challenging: the crystal structures show that the Lewis 
antigens sit in a shallow hydrophobic pocket and that the only 
hydrogen bonds are between the basic residues His85, Argll2, 
and Argl20, and the a 1-2 linked fucose. Few contacts are made 
between LLY^^*^ residues and other carbohydrate components of 
the Le^ and Le^ antigens (Feil etal., 2012). In addition, the struc- 
tural differences between the Le^ and Le^ antigens are minor. 
These characteristics outline the challenges involved in redesign- 
ing LLY^^^ in order to increase its specificity for Le^. Nevertheless, 
the results described herein show that, despite the strong structural 
similarities between Lewis antigens Le^ and Le^ and the shallow 
binding pocket, it is possible to alter the substrate specificity of the 
LLY lectin domain. 

It has previously been demonstrated that LLY^^^wt domain 
binds fucose alone (Farrand etal., 2008; Feil etal., 2012) 
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and that the presence of LLY^^^ enhances the pore-forming 
activity of LLY on platelets. Fucose is a commonly expressed 
carbohydrate, terminally decorating cell surface glycoproteins, 
including CD59 (Rudd etal, 1997; Wheeler etal, 2002), which 
is also the receptor for LLY (Wickham etal, 2011). This sug- 
gests the function of the lectin domain is perhaps as a "capture 
mechanism," slowing down the passing monomers and increas- 
ing the local concentration of LLY in the microenvironment 
of the cell surface. A localized increase in LLY concentration 
would facilitate the oligomerization of LLY into prepores, prior 
to pore formation and consequent cell lysis. In this scenario, 
weak binding affinities for the fucose moieties, as demonstrated 
by SPR, would be beneficial as the monomers would need 
to be readily released for incorporation into assembling pre- 
pore oligomers. Further studies are required to confirm this 
mechanism. 
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